THE WEIGHTED EULER-MACLAURIN FORMULA FOR A SIMPLE 

INTEGRAL POLYTOPE 



JOSE AGAPITO AND JONATHAN WEITSMAN 

Abstract. We give an Euler-Maclaurin formula with remainder for the weighted sum of the values 
of a smooth function on the integral points in a simple integral polytope. Our work generalizes the 
formula obtained in |KSW2| . 



1. Introduction 

The Euler-Maclaurin formula computes the sum of the values of a function / over the integer 
points in an interval in terms of the integral of / over variations of that interval. Khovanskii and 
Pukhlikov |KP11 IKP2j and Kantor and Khovanskii |KK| generalized the classical Euler-Maclaurin 
formula to higher dimensional convex polytopes A which are integral and regular. This formula was 
generalized to simple integral polytopes by Cappell and Shaneson |CS1|ICS2] . and subsequently by 
Guillemin [U] and by Brion-Vergne |BVj . These generalizations involve corrections to Khovanskii's 
formula when the simple polytope is not regular. These formulas are exact formulas, valid for sums 
of exponential or polynomial functions. With the use of the L class associated to the signature 
operator, as in CS1] ICS2j . Karshon, Sternberg and Weitsman gave an Euler-Maclaurin formula 
with remainder for a weighted sum of the values of an arbitrary smooth function on the lattice 
points in a simple integral polytope |KSW1I KSW2 . 

The purpose of this paper is to give a generalization of the Euler-Maclaurin formula with remain- 
der of KSW2 , to allow for more general weightings, including the ordinary, unweighted sum. To 
do this, we use the Hirzebruch formal power series Xq(S) [H| and the weighted polar decomposition 
of[IJ. 

2. Weighted Euler-Maclaurin in one dimension 

Weighted sums in one dimension. Let q be any complex number and f{x) be any function on 
the real line. For any integers a < b, define 

(1) Y?f : = + /(o + 1) + ... + /(&- 1) + ?/(&)• 

[a,b] 

Similarly, for a ray [a, oo) with / having compact support 

(2) Y, 9 f--=qf(a) + f(a + l) + ... 

[a,oo) 

SO 

(3) E'/=£ 5 / -£'-</. 

[a, 6] [a,oo) [b,oo) 

If q = 1 this is the ordinary sum; if q = 1/2 this is the weighted sum of KSW2 . 
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Euler-Maclaurin formulas. The classical Todd function is defined by 



Td(5) 



1 



l_ = 1 - blS+ f2^ 7 s i1 ' 



^ (2n)\ 

n=l v ' 

where b n is the n-th Bernoulli number (we are using the conventions in [Ej). Similarly, the L 
function is given by 



US) 



s / 2 - 1 + f; 4^t.v 2 



tanh(5/2) ^ (2n)l 

Both Td(S') and L(S) are convergent power series for \S\ < 2ir. 



The Hirzebruch function is defined by 



(4) 



1 °° h 

Xq (S) = qTd(S) + (1 - q)Td(-S) = 1 + (q - -)S + £ 7^y S 2n . 

n=l 



It is related to Td and to L by 
(5) Xg (S) = (q- l)S + Td(S) = (q - S + L(5). 

Thus if q = 1 we have xi(S) = Td(5) and if q = 1/2, we get Xl/2(S) = L(S'). 

Let / be a compactly supported function on the real line of class C m where m > 1. The standard 
Euler-Maclaurin formula with remainder for a ray can be written (see for instance |Bj) 



(6) 



y f = Td 2fc (^) r 

[a,ooJ 



f(x)dx 



h=0 



where k = \_m/2\, where Td 2fc (S') denotes the truncation of the power series Td(S') at the 2/c-th 
term, and where the remainder i?^(/) is given by 

B m ({x}) 



poo 

Rm(f) = (-I) 7 "" 1 / P m {x)f {m) {x)dx With P m {x) 
J a 



rn l 



Here, B m (x) is the mth Bernoulli polynomial (see [B]) and {x} = x — \_x\ is the fractional part of 
x. Moreover, the function P m (x) is given by 



(7) 

if m = 2k, and by 

(8) 



P2k{x) 



\fc-l 



E 



2 cos 2imx 



< (27rn) 



2A: 



_ , \ i 2sin27rnx 



if ra = 2fc + l. 

Similarly, denoting the truncation of the power series Xg & t the 2/c-th term by Xq k i we have a 
weighted Euler-Maclaurin formula with remainder for a ray, 



(9) 



E 9 / = xf(|) f°/(^ 

[a,ooJ 



h=0 



where / is a compactly supported function. 
Combining @ and (j2J) we get 



(10) 



r 1.1 1 Ja—hi 

[a,b] 



f{x)dx 



h 1= 



r) f'°° 
dh 2 J b -h 2 h2=0 
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where 

Rt M (f) = R b m(f) ~ RUf) = ("I)™" 1 f P m {x)f (m \x)dx. 

J a 

Given the fact that Xq°(S) is a polynomial whose constant term is 1, we can write this as 

(11) Y J q f=x 2 q k (^-)xfL q (^-)(r f(x)dx- r f(x)dx) +R^ b \f). 

f^J \Ja-h! Jb~h 2 J hl=h2=0 

To make further progress, we note the following symmetry property. 

(12) X f(S)=xf_ q (-S). 

(To see this, observe that Xq {&) is a polynomial with constant coefficients whose constant term 
is 1, whose linear term is (q — 1/2)5 (see ©), and whose other terms are all of even degree 
independent of q.) 

Thus we obtain the following result: 

Theorem 1 (Euler-Maclaurin with remainder for intervals). Let f(x) be a function with m > 1 
continuous derivatives and let k = [m/2\. Then 

fb+h 2 

+ Rt b] (f)- 

h 1 =h 2 =0 



[a,b] 

Note that our argument applies to functions / of compact support. However, for a general 
function / of type C m , the theorem remains true: We need only multiply / by a smooth function 
of compact support which is identically one in a neighborhood of [a,b]. 

If / is a polynomial, Theorem ^ becomes exact when m is greater than the degree of /: 

Corollary 14. Let / be a polynomial. Then 



h 1 =h 2 =o 



(i5) y^x,(^)x,(^) 

[a,b] 



Twisted weighted Euler-Maclaurin for a ray. Consider the "twisted weighted sum" 

oo 

(16) = tf(o) + j>»/(»), 

n>0 rt=l 

where A ^ 1 is a root of unity, say, of order N. 

Let Qm t \(x) be distributions successively defined by 

Qo,a(z) = - J>"5(x-n) 

and 

d f N 
-^Qm,\{x) = Q m ~i,\(x) and J Q m ,\(x)dx = 0. 

These distributions appear in JKS W2 . 
Now, define the polynomial 



N"' A (5) = ( q + — — ) S + Q 2 ,x(0)S 2 + Q 3 ,a(0)5 3 + • • • + Q M (0)S fc , 



A 



1-X / 

for a root of unity A ^ 1. When q = 1/2, we get the polynomial M fc ' A (5) defined in [KSW2] . Since 
N q ' X (S) and M k ' X {S) diff er by (q - l/2)S, adding this term gives the following generalization of 
Proposition 23 in |KSW2j . 
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Proposition 17. Let k > 1 and let f € C 



Then 



(18) 



n>0 17 ' 



f(x)dx 



+ (-1 



,fc-l 



As in |KSW2| . we have the following symmetry property 



(19) 

Remark 20. For A 



N^4 /A (5) 



N 



-S). 



1, we define 



NM(5) 



Q kt \{x)f (k) (x)dx 



xT /2 KS) and Q M = P fc , 



then (HB) boils down to ©. So, (|IBJ| also holds for A = 1. Notice that if A / 1 then N*' A (S) is 
a multiple of S, and that if A = 1 then Ng' A (5) = 1+ a multiple of S. Property (fl9|) continues 
to hold for A = 1 because of the symmetry property (|12jl . Finally, if A = 1, we get the truncation 
Td^S 1 ) of Td(5) for q = 1, and the corresponding truncation at of Td(— 5) for (7 = 0. These 
two expressions differ by S. 



3. Simple polytopes and finite groups associated to them 

In this section we recall various combinatorial and group-theoretic facts about simple polytopes 
which will be needed in our proof of the weighted Euler-Maclaurin formula with remainder. Most 
of this material is taken from and |KSW2j and is included for completeness. 

Let Hi = {x I (ui, x) + [i{ > , Hi £ M} be a half space in W 1 , where u\ € M. n * for i = l,...,d. 
A compact convex polytope A in R n is a compact set which can be written as the intersection of 
finitely many half-spaces 

(21) A = H x n . . . n H d , 

with the smallest possible d, so that the H^s are uniquely determined up to permutation. We order 
them arbitrarily. The facets (codimension one faces) of A are 

Oi = A n dHi , i = l,...,d. 

The vertices of A are all possible codimension n faces obtained by intersections of facets. 

A polytope is called integral if its vertices are in the lattice Z n ; it is called simple if exactly n 
edges emanate from each vertex; it is called regular if, additionally, the edges emanating from each 
vertex lie along lines which are generated by a Z-basis of the lattice Z n . 

For each vertex v of A, let I v C {1, . . . , d} encode the set of facets that contain v, so that 

i G I v if and only if v E <7j. 

The vector itj G M. n * can be thought of as the inward normal to the ith facet of A; a-priori it is 
determined up to multiplication by a positive number. If the polytope A is integral, then the Uj's 
can be chosen to belong to the dual lattice Z n *, and we can fix our choice of the u^s by imposing 
the normalization condition that the u^s be primitive lattice elements, that is, that no u\ can be 
expressed as a multiple of a lattice element by an integer greater than one. 

Assume that A is simple, so that each vertex is the intersection of exactly n facets. For each 
i G I v , there exists a unique edge at v which does not belong to the facet Oi\ choose any vector 
cti tV in the direction of this edge. These vectors form a dual basis to the inward normal vectors n,'s 
and are uniquely determined when the polytope is integral and the u^s are normalized in the sense 
explained above. 
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A "polarizing vector" £ € M. n * is a vector such that (£, aj „) is non-zero for all vertices v and all 
edges i emanating from v. A choice of a polarizing vector £ determine polarized edge vectors oq 
defined by 



(22) 



(23) 



h J Oi.t; if (£, < 0, (unflipped) 
a' = \ 

\—oti,v if > 0, (flipped) 

Let g be any complex number. For each i E I v define 



J. u 

1 - g if af v = -a i; 



l. V ■ 



The tangent cone to A at v is 

(24) C„ = {v + r(x - v) | r>0,ieA} = D + ^ 

Similarly the polarized tangent cone at v is defined by 

(25) Ct = v + J2^>o4,v 

A simple integral orthant C in W 1 is the intersection of n half-planes in general position, 

C = H\ (~l . . . n H n with = {x I {ui,x) + Hi > , fii e R} for i = 1, . . . , n, 

where the Uj's are inward normals to the facets o"j = C n S-ffj of C, which can be chosen to be 
primitive elements of the dual lattice Z n *, and whose vertex v = (~l" =1 {x | (ui, x) + fii = , fii € M} 
is in . This implies that ^ € Z for all 1 < i < n. If ai, . . . , £% is the dual basis to the Uj's, that 
is, 

1 j = i 



J^i, 



(uj,ai) 

then 

(26) C = v + J2~- 

j=l 

We associate a complex number gj to each facet <7j of C, and define the weighting function 

(27) w(x) = n a , 

where ij? denotes the set of facets in C meeting at the face F (the smallest dimensional face in C 
containing x.) If x is in the interior of C, we set w(x) = 1, and if x C, we set w(x) = 0. 
For a simple integral orthant, we consider the weighted sum 

(28) Yl 

cnz™ xecnz" 

where q denotes the n-tuple (qi, ■ ■ ■ q n ) used in the definition of w(x). 

Let q be a complex number. For a polytope, we consider the weighted sum 

(29) E 9 /= E 9 c(a0 /(z), 

Anz n igAnZ" 

where c(x) is the codimension of the smallest dimensional face in A containing x. 
Given q S C, a cone C v and a polarizing vector £ S M ra *, we get the n-tuple q| := (gj v , . . . , qt.,v) 

where q\ v is defined in (|23|) . With this notation, the weighted polar decomposition of JAJ is as 
follows. 
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Theorem 2 (0). For any polarizing vector £, we have 

(30) £ V = £(-i)^ £ q 7, 

where we sum over the vertices v of A and where #v is the number of edge vectors at v that are 
"flipped" by the polarization process h2$fy. 

In obtaining a weighted Euler-Maclaurin formula with remainder for simple integral polytopes, 
we associate certain finite groups to them. We recall some definitions and results from Section §5 
in |KSW2| . 

Let us consider a simple integral orthant C with vertex v € IT 1 , as given in Q26|) . To C we 
associate the finite group 

(31) r :=Z n */Y, Zu i- 

This group is trivial exactly if C is regular. 

Now, let A be a simple integral polytope in R n . For any face F of A, let If denote the set 
of facets of A which meet at F. Because A is simple, the vectors Ui, for i € If, are linearly 
independent. Let Np C M. n * be the subspace 

Np = spanjiij | i G Ip}. 

To each face F of A we associate a finite abelian group Tp- Explicitly, the lattice 

V F =Y^ Zu i <= n f 

is a sublattice of Np nZ n * of finite index, and the finite abelian group associated to the face F is 
the quotient 

(32) r> := (N F nz n *)/v F . 

If F = v is a vertex, this is the same as the finite abelian group associated to the tangent cone C v 
as in (|3Tj) . 

Let E and F be two faces of A with F C E. This inclusion implies that Ip C Ip, and hence 

r B c r>. 

We define a subset T F of Tp by 

(33) t f -.= t f ^ |J r B . 

faces _B such that EZ)F 

Then 

(34) r„= □ r^. 

The map 

(35) A 7J> := e 2^{7,« J ,,) ; for 7 G r fl and j G /wj 
is a well defined character and it is a root of unity. 

Claim 36 ( KSW2 , Claim 61). Ifj€ Tp and j € Ip, then X^j tV is the same for all v 6 F. 
This allows us to define ^,j,F for 7 G Tp and j € If such that 

^7,j,F = A 7 j )U for 7 G and j & Ip, if t> G F. 
Claim 37 ( KS W2j . Claim 62). If 7 GTp and j € I v "^If then Xy* v is equal to one. 
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This allows us to define A 7Jj ^ = 1 when 7 G Fp and j G {1, . . . , d}\Ip. Then 

(38) ^y,j,F = ^7,j> for 7 £ Tj? and 1 < j < d, if v G F 
and 

(39) A 7ji F = 1 for 7 G r> if j i>. 
Claim 40 C jKSW2j . Claim 65). 1/76^ and j G I F , then A 7J - F ^ 1. 

4. Weighted Euler-Maclaurin with remainder for simple integral polytopes. 

Let r be the finite group (|31[) associated to a simple integral orthant C (as defined in (|26j0 with 
vertex at d £ Z". The map 7 1— > e^ 7 ' 3 ^ defines a character on T f |KSW2| . Lemma 52) whenever 
x G an d this character is trivial if and only if x G Z ra . 

By a theorem of Frobenius, 

1 ^jMi,*)-! 1 if*e^ n 



if x Z n 



for all x G ^]Zaj. Then, for any vector q = (qi, . . . ,q n ) G C n and any function f(x) compactly 
supported on W 1 , 

cnz™ x y ' 7 er / 

(41) = ^EE qe2m<7 '^) 

where we sum over all 

(42) x = v + miai + . . . + m n a n , 

with the rrij's being non-negative integers. 

The simple integral orthant C is the image of the standard orthant O = YYt=\ ^>0 i n ^ n under 
the affine map 

(ti, . . . ,t n ) H-> V + E^ a *- 

This map sends the lattice Z n onto the lattice Y^^ a j- Let us concentrate on one element 7 G T. 
Because v G Z n , from (fl2*|) we get 

n 

e 2^< 7 ,*) = JJ A ^ where x . = e 2ni( 7 , aj ) ^ 

so that the inner sum in H41[) becomes 

(43) YT e2Ml ' X) f( x ) = E 91 A r • • • E "^n" 9(mx,. . .,m n ), 

x mi>0 m n >0 

where 

5(7711, . . . , m n ) = f(v + mi«i + . . . + m n a n ). 
Given q G C and > 1, we had the twisted remainder formula (see (J18|0 

Q /-OO />0O 

E A m 5 (m) = N/' A (— ) / g(t)dt + / Q k ,x(t)g {k) (t)dt 

for all compactly supported functions g(x) of type C fc , where k > 1, where A is a root of unity, and 

A: A 

where iV g ' is a polynomial of degree < k. 
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Iterating this formula, the sum in (|43|) can be written as 

d 



d f 

n N t' Al (^) y g(h,...,t n )dt 1 ---dt n + R% k (\ 1 ,...,\ n ;g), 



i_1 O(h) 
where 

0(/ii, . .. ,/i n ) = {(*!, ... ,i n ) | ij > -/tj for all i} 
denotes the "dilated" standard orthant, and where the remainder is given by 

(44) R^ k (Xi,...,X n ;g):= £ (-l)(*-i)(»HD 

/C{l,...,n} 



of U 



ie/ ' 0(h) Hi Hi 

with 

g(h, ...,t n ) = f(v + t x ai + . . . + t n a„). 
Performing the change of variable given by the transformation 

L: (ti, . . . ,t n ) i-> x = v + tiai + . . . + t n a n , 

whose Jacobian is l/|r|, and substituting back into l|41|l. we get 



h=0 



(45) 



E q / = EII N ^(^) / /(*)* 

cnz« 7 eri=i 1 c J {h) 



h=0 

where C(h\, . . . , /i n ) denotes the image of the "dilated" standard orthant 0(hi, . . . , h n ) under the 
affine transformation L, and where the remainder is given by 

(46) R° k (f) := ^ E <fc( A 7,i> ■ ■ ■ > A 7 ,„; LV), 

where 

A 7j := e 2 ™<^>. 



Let now A be a simple integral polytope, given by 1)21(1. Choose a polarizing vector for A and 
let Cf, denote the polarized tangent cones. We can consider the "dilated polytope" A(/ii, . . . , h^), 
which is obtained by shifting the ith facet outward by a "distance" hi. More precisely, 

d 

A(/i) = I (ui,x) + ^ + hi > 0} where h = (hi, . . . , h d ). 

i=i 

Then A(h) is simple if h is sufficiently small. The polar decomposition of A(/i) involves "dilated 
orthants" . However, dilating the facets of A outward results in dilating some facets of Cf, inward 
and some outward. Explicitly, for i € I v = {ii, . . . ,i n }, the inward normal vector to the ith facet 
of C| is 



(47) u 



l,V 



•r |J _ 

Hi II OL^ v — (%i,v 
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Hence, the dilated orthants that occur on the right hand side of the polar decomposition of A(/i) 
are C|(/if , . . . , h\ n ), where 



(48) 



h 



hi if a^ v = a i)V 



-hi if a] 



Because the inward normals to the facets of C| are given by (|47l) . the dual basis to these vectors 
is a| v , i G I v , and the roots of unity that appear in the Euler-Maclaurin formula for C{, are then 



A" ■ = e 



2ni 



(^«L) = J A W> ' lia L = a i,v 



XZL, if a 



(49) 

Let k > 1 be an integer. For any compactly supported function / on R n of type C nk , we then 
get from (EOJ) and (gSJ) 



Anz n f 



11 ■■' 



iiv dh* ' 



f(x)dx 



CUh s ,...M ) 



h=0 



+ <*(/), 



where the remainder is given by 

(si) RfM) ■■= Ec- 1 )*^ *(/) 



and where the fife's are given in (|15|) . Note that either h\ v = hi, Xi /iv = Xj,i, v and q\ v = q, or 



-IH, A 1} i }V 



X \ v and q?^ = 1 — q. By the symmetry property (|T9|) . we have 



(52) 



4, q { dh~ 



Now, for j g" /„, because A 7Jil) = 1 (see (jSHJ)), we have N g fe ' A ™>(g|-) = 1+powers of ^f-. Also 

for j g" I v , the cone Cj(/tj- /if n „) is independent of /ij. These facts together with Formula (|52|) 

imply that (|ST)|) is equal to 



(53) 



7er„ i=i 



/i=0 



Because Xjj t F = A 7Jil) whenever v £ F (see (|3*5|) ). we can define 

d 

(54) 

3=1 

and we have 



(55) N g, 7 ,F = N g, 7 ,i> whenever 7 € r F and v G F, 

where we identify 7 G with its image under the inclusion map <— > T r 
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Then (|53~)) is equal to 



(56) 



B- 1 )*" E K„ I 



f(x)dx 



+ <*(/) 



ft=0 



'C»(ft». .....ft* ) 



f(x)dx 



+ R£k(f)> 



h=0 



by (|34|) and (|55j). In the interior summation we may now add similar summands that correspond 
to v F. These summands make a zero contribution to (|56j) for the following reason. If v F 
then there exists i £ If^I v . Because i £ I v , the cone Cj^j^ • ■ ■ fojj n „) is independent of /ij. So 
it is enough to show that F is a multiple of -M-. But because 7 G T F and i £ If, we have 



A 7) i,ir ^ 1. (See Claim E01) By Remark HOI this implies that N,' 7, ' ,F (^-), which is one of the 
factors in N^ 7 ^, is a multiple of Hence, (|56|) is equal to 



(57) EE N ^E(-W 

F 7gr t all?; 17 C 



iciM ,...,ft» ) 



f{x)dx 



h=0 



E E N ",7,F / f(x)dx 
F 7Gr ^ ^(ft) 



R £k(f) 



R tkU)- 



h=0 



We have therefore proved our main result: 

Theorem 3. Let A be a simple integral polytope in W 2 . Let k > 1 and let f € C™ fc (lR n ) be a 
compactly supported function on W l . Choose a polarizing vector for A. Then 



E^ = E E n Lf f f{*)dx 

AnZ" f 7g r b JA(ft) 



+ <*(/) 



ft=0 



where ^ are differential operators defined in (|54|1 and where the remainder Rg^if) ^ s 5* wen 
equation \51)) . Moreover, the differential operators ^ are 0/ order < k in each of the variables 
hi, . . . ,hd- Also, the remainder can be expressed as a sum of integrals over orthants of bounded 
periodic functions times various partial derivatives of f of order no less than k and no more than 
kn. This remainder is independent of the choice of polarization and is a distribution supported on 
the polytope A. 

In particular, if q = 1/2, we get Theorem 2 of |KSW2] . while if q = 1, we have a formula for the 
ordinary, unweighted sum. Notice that if q = 0, we also have an unweighted sum but only over the 
interior lattice points in the polytope. 



5. Estimates on the remainder and an exact Euler-Maclaurin formula for 

polynomials 

In order to derive a formula for polynomials from Theorem |HJ we first require an estimate on 
the remainder term R^ k (f). Returning to the definition of functions Q m ,\{x) in KSW2 j; we see 
that Qm,x(x) is a periodic function on R. It follows that Q m ,x{x) is bounded. Since the operators 
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Ng' A (^) are differential operators of order k, the definition of the remainder R^ k (Xi, • • • , A n ; /) 
(see (144(1 ) shows that 

(58) l< fc (Ai, • • • , A n ; /)| < JT q • sup i;; . . . dtf\ Ll(Rn) , 

where the supremum is taken over all n-tuples ■ ■ ■ ,j n } with k < j\ + ■ • • + j n < nk. 

The definition of R^ k and equation (|58|l then give the same estimate for the remainder on the 
polytope, of course with a different constant. 

Proposition 59. The remainder term in Theorem^ can be estimated by 

\R£ k (f)\ < K q (k, A) • SUp; ; . ; , ;/>[ . . . ^/| il(K n), 

where the supremum is taken over all n-tuples {ji, ■ ■ ■ ,j n } with k < j\ + ■ ■ ■ + j n < nk. 

The estimate in Proposition 1591 implies 
Proposition 60. Let p be a polynomial on M n , and choose k > deg p + n + 1. Then 



\C\7S 1 F „,^-rb J A(h) 



f h=0 

Corollary 61. Let p be a polynomial and suppose that the polytope A is regular. Then 

d ^ 

N Ix 

1 ' nil , i 



q p = tl^(^h) [ p^ da 
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